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ABSTRACT 

General relativistic spectra from accretion disks around rotating neutron stars in the appropriate 
space-time geometry for several different equation of state, spin rates and mass of the compact object 
have been computed. The analysis involves the computation of the relativistically corrected radial 
temperature profiles and the effect of Doppler and gravitational red-shifts on the spectra. Light bending 
effects have been omitted for simplicity. The relativistic spectrum is compared with the Newtonian one 
and it is shown that the difference between the two is primarily due to the different radial temperature 
profile for the relativistic and Newtonian disk solutions. 

To facilitate direct comparison with observations, a simple empirical function has been presented 
which describes the numerically computed relativistic spectra well. This empirical function (which has 
three parameters including normalization) also describes the Newtonian spectrum adequately. Thus the 
function can in principle be used to distinguish between the two. In particular, the best fit value of 
one of the parameters (/3-parameter) w 0.4 for the Newtonian case, while it ranges from 0.1 to 0.35 
for relativistic case depending upon the inclination angle, EOS, spin rate and mass of the neutron star. 
Constraining this parameter by fits to future observational data of X-ray binaries will indicate the effect 
of strong gravity in the observed spectrum. 

Subject headings: X-rays:binaries-X-rays:spectra -stars:neutron -stars:rotation -relativity:general 



1. INTRODUCTION 

X-ray binaries are believed to harbor black holes or 
weakly magnetized neutron stars with an accretion disk. 
The X-ray emission arises from the hot (« 10 7 K ) in- 
nermost region of the disk. In the case of a neutron star 
there will emission, in addition, from a boundary layer be- 
tween the accretion disk and neutron star surface. Since 
the observed emission arises from regions close to compact 
object, these sources are possible candidates for studying 
strong field gravity. 

In the standard theory (Shakura & Sunyaev 1973), the 
accretion disk is assumed to be an optically thick Newto- 
nian one. In this model, the local emergent flux (assumed 
to be a blackbody) is equated to the energy dissipation 
at a particular radial point in the disk. The observed 
spectrum is then a sum of black body components aris- 
ing from different radial positions in the disk. General 
relativistic effects modify this Newtonian spectrum in two 
separate ways. Firstly, the local energy dissipation at a 
radial point is different from the Newtonian disk, giving 
rise to a modified temperature profile. Secondly the ob- 
served spectrum is no longer a sum of local spectra because 
of effects like Doppler broadening, gravitational red-shifts 
and light bending. Modified spectra, incorporating these 



effects, but with different approximations have been com- 
puted by several authors (e.g. Novikov & Thorne 1973; 
Asaoka 1989) for accretion disks around rotating (Kerr) 
black holes. These computations confirm the expected re- 
sult, that the relativistic spectral shape differs from the 
Newtonian one by around 10%. Thus, for comparison with 
observed data with systematic and statistical errors larger 
than 10%, the Newtonian approximation is adequate. Ebi- 
sawa, Mitsuda and Hanawa (1991) showed that for typical 
data from the Ginga satellite, the relativistic spectrum 
cannot be differentiated from the Newtonian disk spec- 
trum. They also found that the relativistic spectrum is 
similar in shape (at the sensitivity level of Ginga ) to the 
Comptonized model spectrum. Although, Ginga was not 
sensitive enough to distinguish between the different spec- 
tra, better estimates of fit parameters like accretion rate 
and mass of the compact object were obtained when the 
data was compared to relativistic spectra rather than the 
standard Newtonian one. 

The present and next generation of satellites (e.g. 
ASCA, RXTE, Chandra, XMM, Constellation-X) with 
their higher sensitivity and/or larger effective area than 
Ginga are expected to differentiate between relativistic 
and Newtonian spectra from low mass X-ray binaries 
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(LMXB) and black hole systems. However, as pointed 
out by Ebisawa, Mitsuda and Hanawa (1991), the pres- 
ence of additional components (e.g. boundary layer emis- 
sion from the neutron star surface) and smearing effects 
due to Comptonization may make the detection ambigu- 
ous. Nevertheless, the detection of strong gravity effects 
on the spectra from these sources, will be limited by the 
accuracy of theoretical modeling of accretion disk spec- 
tra rather than limitations on the quality of the observed 
data. Thus, it is timely to develop accurate relativisti- 
cally corrected spectra for comparison with present and 
future observations. Apart from the importance of detect- 
ing strong gravity effects in the spectra of these sources, 
such an analysis may also shed light on the geometry and 
dynamics of innermost regions of accretion disks. 

Novikov & Thorne (1973) and Page & Thorne (1974) 
computed the spectra of accretion disks around rotating 
(Kerr) black holes. This formalism when directly applied 
to rotating neutron stars provide only a first order esti- 
mate: the absence of an internal solution in the case of 
Kerr geometry, makes it difficult to obtain, in a straight- 
forward manner, the coupling between the mass and the 
angular momentum of the central accretor. This coupling 
depends on one hand on the equation of state of neu- 
tron star matter and on the other on the proper treat- 
ment of rotation within general relativity. Equilibrium 
configurations of rapidly rotating neutron stars for real- 
istic equations of state have been computed by several 
authors (Bonazzola & Schneider 1974; Friedman, Ipser & 
Parker 1986; Cook, Shapiro & Teukolsky 1994; Stergioulas 
& Friedman 1994; Salgado et aL1994a,b; Datta, Tham- 
pan & Bombaci 1998). One crucial feature in all these 
calculations is that the space-time geometry is obtained 
by numerically and self-consistently solving the Einstein 
equations and the equations for hydrodynamic equilibrium 
for a general axisymmetric metric. With the aim of mod- 
eling spectra of LMXBs, we attempt, in this paper, to 
compute the spectrum of accretion disks around rotating 
neutron stars within such a space-time geometry. This is 
particularly important since LMXBs are old (population 
I) systems and the central accretor in these systems are ex- 
pected to have large rotation rates (Bhattacharya & van 
den Hcuvcl 1993 and references therein) 

Computation of the spectra is numerically time consum- 
ing and hence direct fitting to the observed data is imprac- 
tical. For the sake of ease in modeling, we also present in 
this paper, a simple empirical analytical expression that 
describes the numerically computed spectra. As shown 
later, the same expression (which has three parameters 
including normalization) can also describe the Newtonian 
spectra. In particular, the value of one of the parameters 
(called (3 parameter here) determines whether the spec- 
trum is relativistically corrected or not. This will facil- 
itate comparison with observational data since only this 
/3-parameter has to be constrained to indicate the effect of 
strong gravity in the observed spectrum. 

The next section describes the method used to compute 
the spectra. In §3 the results of the computation and the 
empirical fits are shown. §4 is devoted to discussion and 
summary. 

2. SPECTRAL COMPUTATION 



The disk spectrum is expressed as: 



F(E oh ) = (l/E oh ) J I oh (E oh )dn oh 



(1) 



where the subscript 'ob' denotes the quantity in observer's 
frame, the flux F is expressed in photons/sec/cm 2 /keV, E 
is photon energy in keV, I is specific intensity and Q is the 
solid angle subtended by the source at the observer. 

As I/E 3 remains unchanged along the path of a photon 
(see for e.g., Misner et al. 1973), one can calculate J b, 
if J em is known (hereafter, the subscript 'em' denotes the 
quantity in emitter's frame). We assume the disk to emit 
like a diluted blackbody, so I cm is given by 

I cm = (l/f 4 )B(E cmi T c ) (2) 

where / is the color factor of the disk assumed to be inde- 
pendent of radius (e.g. Shimura & Takahara 1995). B is 
the Planck function and T c (the temperature in the central 
plane of the disk), is related to the effective temperature 
T e g through the relation T c = /T e g. The effective temper- 
ature, T e g is a function of the radial coordinate r and for 
a rotating accretor, is given by (Page & Thorne 1974) 

T cff = (F/a) 1 / 4 

where 



(3) 
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tl K l)l, r dr (4) 



is the flux of energy from the disk in orbiting particle's 
frame. The factor a is the Stcphan-Boltzmann constant, 
r m is the disk inner edge radius, E, I are the specific en- 
ergy and specific angular momentum of a test particle in a 
Keplerian orbit and is the Keplerian angular velocity 
at radial distance r. In our notation, a comma followed by 
a variable as subscript to a quantity, represents a deriva- 
tive of the quantity with respect to the variable. Also, in 
this paper, we use the geometric units c = G = 1. 

The quantities E oa and E cm are related through the ex- 
pression E cm = E Q b(l+z), where (1+z) contains the effects 
of both gravitational redshift and Doppler shift. For a gen- 
eral axisymmetric metric (representing the space-time ge- 
ometry around a rotating neutron star), the factor (1 + z) 
is expressed as (see for e.g. Luminet 1979) 

1 + z = (1 + Q, K b sin a sin i)(-g tt - 2Vt K g t4> - Vt^g^) \&) 

where the g^s are the metric coefficients and t and <f> 
are the time and azimuthal coordinates. In the above ex- 
pression (which includes light-bending effects), i is the 
inclination angle of the source, b the impact parameter of 
the photon relative to the line joining the source and the 
observer and a the polar angle of the position of the pho- 
ton on the observer's photographic plate. For the sake of 
illustration and simplicity in calculations, we neglect light 
bending. We thus write b sin a = r sin <f> and 



rdrd(j) cos i 



(6) 



where D is the distance of the source from the observer. 

The space-time geometry around a rotating neutron 
star can be described by a general axisymmetric, station- 
ary metric (see for e.g. Komatsu, Eriguchi & Hachisu 
1989). Assuming the matter to be a perfect fluid, and the 
metric to be asymptotically flat, the Einstein field equa- 
tions reduce to three non-homogeneous, second-order, 
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coupled differential equations and one ordinary differ- 
ential equation in terms of the energy density and the 
pressure of the high-density neutron star matter (Cook, 
Shapiro & Teukolsky 1994). An important input in solv- 
ing these equations is the equation of state (EOS) of high 
density matter composing the neutron star. Assuming 
rigid rotation, we solve these equations numerically and 
self-consistently for four representative EOS models: (A) 
Pandharipande (1979) (hyperons), (B) Baldo, Bombaci & 
Burgio (1997) (AV14 + 3bf), (C) Walecka (1974) and (D) 
Sahu, Basu & Datta (1993). The values of stiffness param- 
eters for each of these EOS models are widely different and 
increases from (A) to (D). We, therefore, expect the results 
of our computations to be of sufficient generality. 

The solution of the field equations yield the metric co- 
efficients, numerically, as functions of r and 8. Using 
these metric coefficients, it is straightforward to calcu- 
late the structure parameters of rapidly rotating neutron 
stars. The details of these calculations are given in Datta, 
Thampan & Bombaci (1998) and references therein. The 
quantities: r- mi E, I and f2x are obtained by solving the 
equation of motion of material particles within the space- 
time geometry given by the above metric (Thampan & 
Datta 1998; Bhattacharyya et aL2000). For our purpose 
here, we compute constant gravitational mass sequences 
whose rotation rates vary from zero to the centrifugal mass 
shed limit (where gravitational forces balance centrifugal 
forces). For realistic neutron stars, the inner radius r; n 
may be located either at the marginally stable orbit or 
the surface of the neutron star depending on its central 
density and rotation rate (Thampan & Datta 1998; Bhat- 
tacharyya et al.2000), having important implications for 
the gravitational energy release as well as the temperature 
profiles of accretion disks. For the procedure of calculat- 
ing T e fF, for rapidly rotating neutron stars, considering the 
full effect of general relativity, we refer to Bhattacharyya et 
al., 2000. These authors have also shown (in their Fig. 2) 
that the difference between Newtonian temperature profile 
and general relativistic temperature profile is substantial 
at the inner portion of the disk. As will be shown herein, 
it turns out that this is the major reason for the difference 
between Newtonian and general relativistic spectra at high 
energies. 

To summarize this section, we calculate the accretion 
disk spectrum using equation 1, taking the radial integra- 
tion limits as ri n and r ou t and the azimuthal integration 
limits as and 2n. We choose a very large value (« 10 5 
Schwarzschild radius) for r out . 

3. RESULTS 

To illustrate the differences between the relativistic and 
Newtonian spectra, we show in Figure 1, the computed 
relativistic spectrum (solid line) and the Newtonian spec- 
trum (dotted line) for the same parameters. The Newto- 
nian spectrum is the spectrum expected from a standard 
non-relativistic disk (Shakura & Sunyaev 1973) but with 
the intensity and the effective temperature modified by the 
color factor (equations 2 and 3). In order to isolate the 
different contributions, we have also plotted in Figure 1, 
the theoretical spectrum arising from relativistic tempera- 
ture profile but without the effect of Doppler/gravitational 
red-shift (dashed line). The relativistic spectrum is under- 



luminous compared to the Newtonian one at high energies 
- this is primarily because of the the difference in the ra- 
dial temperature profile (Bhattacharyya et al. 2000). The 
difference between the two spectra is nearly 50% at 2 keV. 
We emphasize here, that such high difference is only true 
when both the spectra are calculated for the same disk 
parameters. If the Newtonian spectra is calculated for 
slightly different values of disk parameters (e.g. accretion 
rate, inclination angle, distance to the source) the average 
discrepancy between the two spectra will be less (Ebisawa, 
Mitsuda and Hanawa 1991). 

In order to facilitate comparison with observations, we 
introduce a simple analytical expression which empirically 
describes the computed relativistic (and Newtonian) spec- 
tra. 

S { (E) = S E-^(^re X p(-^) (7) 

where, 7 = — (2/3)(l + Ep/E a ), E a , /3 and S are parame- 
ters and E is the energy of the photons in keV. Sf(E) is in 
units of photons/sec/cm 2 /keV. To compare this empirical 
function with the computed spectra, we use a reduced x 2 
technique. In particular, we define a function 

2 _ J_ sr^ ( Sc( E i) - Sf(Ej) 2 

N ii °- ls c{Ei) ' ' 

where S C (E) is the computed spectra. The spectra are 
divided into N logarithmic energy bins. We have chosen 
the range of energy used in calculating x 2 to be dependent 
on the location of the maximum of the energy spectrum 
( ES C (E)) which is typically at 2 keV. The minimum en- 
ergy is set to be one hundredth of this value (typically 0.02 
keV) while the maximum is set at ten times (typically 20 
keV). x 2 is fairly insensitive to the number of energy bins; 
we take N — 200. For each S C (E) the best fit parameters 
(E a , P and S ) are obtained by minimizing x 2 ■ Figure 2 
shows the relativistic spectra for three different inclination 
angles (solid lines) and the corresponding empirical fits us- 
ing equation (1) (dotted lines). The minimum x 2 obtained 
while fitting these spectra was < 0.1, which means that the 
average discrepancy is less than 3 %. This is also true for 
other disk parameters and EOS considered in this work. 
Thus the empirical function (eq. [1]) is a reasonable ap- 
proximation to the computed relativistic spectra. It also 
describes the Newtonian spectra to a similar degree of ac- 
curacy. 

The S parameter in equation (1), is the normalization 
factor and is independent of the relativistic effects. It de- 
pends only on the mass of the star (M), accretion rate 
(M), distance to the source (D), color factor (f) and in- 
clination angle (i) i.e S Q oc M 2 / 3 f- 4 / 3 M 1 / 2 D- 2 cos i. The 
E a parameter (which is in units of keV) describes the high 
energy cutoff of the spectrum. Its dependence on the space 
time metric and inclination angle is complicated but it 
scales as E a oc M 1 / 4 /. The /3-parameter depends only on 
the space time metric and not on accretion rate, distance 
to the source or color factor. This makes the /3-parameter 
useful as a probe into the underlying space-time metric. 
We show in Figure 3, the variation of minimum x 2 (i- e - 
minimized w.r.t. to parameters E a and S only) as a 
function of the /3-parameter for the three spectra shown 
in Figure 2 and for the Newtonian one. For the Newto- 
nian case the minimum x 2 occurs for (3 w 0.4 while it 
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is lower for the relativistic cases. For example, consider 
the relativistic spectrum for parameters listed in Figure 1 
and for i = 30° (line marked as 3 in Figure 3). If this 
spectrum is fitted with the empirical function the mini- 
mum x 2 = 0.05 (corresponding to an average discrepancy 
of 2%) and the best fit /3-parameter is (3 ~ 0.25. For 
a Newtonian /3-parameter value of « 0.4, the minimum 
X 2 increases to 0.1, corresponding to an average discrep- 
ancy of more than 3%. Thus the empirical function can 
resolve the difference between the Newtonian and the rel- 
ativistic one at the 10% level. For an observed spectrum 
fitted using the empirical function, if the best fit range 
of /3-parameter excludes the Newtonian value of 0.4, that 
would strongly indicate that the spectrum has been modi- 
fied by strong gravitational effects. Since the /3-parameter 
deviation from 0.4 increases with inclination angle, nearly 
edge on disks are more promising candidates for detecting 
the presence of strong gravity. To show the robustness of 
this result we show in Figure 4, 5 and 6 the variation of 
the best fit /3-parameter with i for different EOS, mass of 
the central object and spin rates respectively. For all these 
cases the best fit /3-parameter is less than 0.4. Parameter 
E a is useful to determine the accretion rate. However, it 
also depends on the metric and inclination angle. We show 
this dependence in Figure 7. 

4. SUMMARY AND DISCUSSION 

In this paper, we have computed relativistic spectra 
from accretion disks around rotating neutron stars for the 
appropriate space-time geometry. Several different EOS, 
spin rates and mass of the compact object have been con- 
sidered. The Doppler and gravitational effects on the spec- 
tra have been taken into account, while light bending ef- 
fects have been omitted for simplicity. The spectrum dif- 
fers from the Newtonian one, with the main difference be- 
ing due to the different radial temperature profile for the 
relativistic and Newtonian disk solutions. 

A simple empirical function has been presented which 
describes the numerically computed relativistic spectra 
well. This will facilitate direct comparison with obser- 
vations. The empirical function (eq. [7]) has three param- 
eters including normalization. Another important advan- 
tage of this function is that it also describes the Newtonian 
spectrum adequately, and the value of one of the parame- 
ters (/3-parameter) distinguishes between the two. In par- 



ticular, the best fit /3-parameter ss 0.4 for the Newtonian 
case, while it ranges from 0.1 to 0.35 for relativistic case 
depending upon the inclination angle, EOS, spin rate and 
mass of the neutron stars. 

In principle, for sufficiently high quality data, the effects 
of strong gravity on the disk spectrum can be detected 
using this empirical function as a fitting routine and con- 
straining the 6-paramcter. However, it must be empha- 
sized that there are several reasons why this may not be 
possible. There could be systems which have additional 
components in the X-ray spectra; for example boundary 
layer emission from the neutron star surface. Uncertain- 
ties in modeling these additional components may lead to 
a wider range in the best fit /3-parameter. Thus accu- 
rate spectra of the boundary layer (with relativistic cor- 
rections) is also needed for modeling these systems . More- 
over, X-rays could be emitted from hotter regions (e.g an 
innermost hot disk or a corona) giving rise to a Comp- 
tonized spectra instead of the sum of local emission as- 
sumed here. In this case, the empirical fit will probably not 
describe the observed data well. It has been assumed here 
that the color factor is independent of radius. Shimura 
and Takahara (1995) have shown from numerical compu- 
tation that this could be the case for an accretion disk in 
a Schwarzschild metric. Apart from the fact that this was 
done for Schwarzschild metric, their numerical computa- 
tion depends on the vertical structure of the disk which 
in turn depends on the unknown viscosity mechanism in 
the disk. If the color factor has a radial dependence, the 
spectral shape might change which may be confused to be 
a relativistic effect. 

Despite these caveats the method described in this paper 
will be a step forward in the detection of strong gravity 
effects in the spectra of X-ray binaries. Future compar- 
ison with high quality observational data, will highlight 
the theoretical requirements that have to be met, before 
concrete evidence for strong gravity are detected in these 
systems and the enigmatic region around compact objects 
is probed. 
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Pijush Bhattacharjec for encouragement. 
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Fig. 1. — General relativistic spectrum (solid line) for a neutron star configuration with mass M = IAMq, spin rate Q* = 0, distance to 
the source D = 5 kpc, inclination angle i = 30°, accretion rate M = 10 18 g/sec and color factor / = 2. Dashed line: the spectrum expected 
from a source with the same disk parameters but without the relativistic effects (Newtonian spectrum). Dotted line: The spectrum for the 
same disk parameters but without the effect of Doppler and gravitational red-shifts (i.e. z is set to zero). The EOS model (B) is used here. 
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Fig. 2. — The relativistic spectra for three different inclination angles ( i = 0°, 30°, 60° ) with rest of the parameters same as in Figure 1 (solid 
lines). Dashed lines: empirical fit to the relativistic spectra using equation (7). The minimum \ 2 = 0.073,0.049,0.026 for i = 0°,30°,60° 
respectively. 
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Fig. 3. — Variation of minimum \ (i- e - minimized w.r.t to parameters E a and So) with parameter /3. Curves marked 2, 3 and 4 correspond 
to the spectra shown in Figure 3 for i = 0°, 30°, 60° respectively. Curve marked 1 is for the Newtonian spectra shown in Figure 1. 
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Fig. 6. — Variation of the best-fit /3-paramcter with inclination angle for different neutron star spin rates. Curve 1: Q t = radians/sec, 
curve 2: f2* = 2044 radians/sec, curve 3: Q* = 7001 radians/sec (mass-shed limit). The values of the other parameters are as in Figure 1. 
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